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ABSTRACT 

We perform a detailed study on the dynamics of a relativistic blast wave with 
the presence of a long-lived reverse shock (RS). Although a short-lived RS has 
been widely considered, the RS is believed to be long-lived as a consequence of 
a stratification expected on the ejecta Lorentz factors. The existence of a long- 
lived RS makes the forward shock (FS) dynamics to deviate from a self-similar 
Blandford-McKee solution. Employing the "mechanical model" that correctly 
incorporates the energy conservation for such blast waves with a long-lived RS, 
we present an accurate solution for both the FS and RS dynamics. We conduct 
a sophisticated calculation of the afterglow emission. Adopting a Lagrangian 
description of the blast wave, we keep track of an adiabatic evolution of numer- 
ous shells between the FS and RS. An evolution of the electron spectrum is also 
followed individually for every shell. We then find the FS and RS light curves by 
integrating over the entire FS and RS shocked regions, respectively. In particular, 
we make use of an analytic expression for observed spectral flux, which we derive 
here in terms of an observed frequency and observer time. Exploring a total 
of 20 different ejecta stratifications, we explain in detail how a stratified ejecta 
affects its blast wave dynamics and afterglow light curves. We show that, while 
the FS light curves are not sensitive to the ejecta stratifications, the RS light 
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curves exhibit much richer features, including steep declines, plateaus, bumps, 
re-brightenings, and a variety of temporal decay indices. These distinctive RS 
features may be observable if the RS has higher values of the micophysics parame- 
ters than the FS. We discuss possible applications of our results in understanding 
the GRB afterglow data. 

Subject headings: gamma-ray burst: general - - radiation mechanisms: non- 
thermal — shock waves 



1. Introduction 

The central engine of a gamma-ray burst (GRB) ejects a relativistic outflow (called 
an ejecta) with high Lorentz factors. As the ejecta interacts with a surrounding ambient 
medium, a relativistic blast wave develops. The blast wave consists of two shock waves: the 
forward shock (FS) wave sweeping up the ambient medium and the reverse shock (RS) wave 
propagating through the ejecta. The shocked ambient medium is separated from the shocked 
ejecta by a contact discontinuity (CD), and a compressed hot gas between the FS and RS is 
called a "blast". 

Without an extended activity of the central engine, the RS is expected to be short- 
lived if the ejecta is assumed to have a constant Lorentz factor r e j. The RS vanishes as 
it crosses the end of the ejecta. The blast wave then enters a self-similar stage where the 
FS dynamics is described by the solution of Blandford & McKee (1976) (hereafter BM76). 
This FS emission has been believed to be the main source of the long-lasting, broad-band 
afterglows (Meszaros & Rees 1997; Sari et al. 1998). The short-lived RS emission would be 
then important only briefly in the early afterglow phase. Thus, it was proposed to explain 
a brief optical flash detected in some GRBs (Meszaros & Rees 1997,1999; Sari & Piran 
1999a,b). The dynamical evolution of such a short-lived RS with a constant r ej was studied 
analytically (Sari & Piran 1995; Kobayashi 2000), under the assumption of an equality of 
pressure across the blast wave. 

However, a general view on the structure of the ejecta should include the possibility 
that the ejecta emerges with a range of the Lorentz factors. The shells with lower Lorentz 
factors gradually "catch up" with the blast wave as it decelerates to a comparable Lorentz 
factor. Therefore, the RS wave is believed to be long-lived in general. An example with a 
long-lived RS, where a power-law ejecta interacts with a power-law ambient medium, was 
studied analytically by assuming a constant ratio of two pressures at the FS and RS (Rees 
& Meszaros 1998; Sari & Meszaros 2000). 
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The structure or stratification of the ejecta and ambient medium could be in fact even 
more general. There is no reason why it should take only a constant or power-law profile. 
Uhm (2011) (hereafter Ull) presented a semi-analytic formulation for this class of general 
problems where the ejecta and ambient medium can have an arbitrary radial stratification. 
Ull takes into account a radial spread-out and spherical expansion of such a stratified ejecta 
and finds which shell of this evolved ejecta gets passed by the RS at a certain time and 
radius. Ull then finds the dynamics of the blast wave with a long-lived RS, by employing two 
different methods: (1) an equality of pressure across the blast wave (mentioned above) and 
(2) the "mechanical model" (Beloborodov & Uhm 2006). Ull shows that the two methods 
yield significantly different dynamical evolutions and demonstrates that the method (1) does 
not satisfy the energy conservation law for an adiabatic blast wave while the method (2) 
does. The mechanical model does not assume either an equality of pressure across the blast 
wave or a constant ratio of two pressures at the FS and RS. It shows that the ratio of two 
pressures should in fact evolve in time as the blast wave propagates. 

Besides these theoretical considerations, recent early afterglow observations led by Swift 
revealed a perplexing picture regarding the origin of GRB afterglows. In contrast of a simple 
power-law decay feature as expected from the standard afterglow theory, the X-ray data 
show more complicated features including initial rapid declines, plateaus, and flares (e.g. 
Tagaliaferri et al. 2005; Burrows et al. 2005; Nousek et al. 2006; O'Brien et al. 2006; 
Chincarini et al. 2007) that reveal rich physics in the early afterglow phase (Zhang et al. 
2006; Zhang 2007). More puzzlingly, some GRBs show clear chromatic behaviors of the 
X-ray and optical afterglows (e.g. Panaitescu et al. 2006; Liang et al. 2007, 2008). It is now 
evident that the FS alone cannot interpret the broad-band afterglow data for the majority 
of GRBs. 

Uhm & Beloborodov (2007) and Genet et al. (2007) independently showed that the 
RS-dominated afterglow flux could reproduce some observed afterglow features, given the 
assumption that the FS emission is suppressed. In this paper, we study in great detail 
the dynamics and afterglow light curves of GRB blast waves with a long-lived RS. The 
purpose is to investigate how different ejecta stratifications affect their blast wave dynamics 
and afterglow light curves. We explore various types of the ejecta stratification and unveil 
that there exists a whole new class of the blast wave dynamics with a rapid and strong 
evolution of the RS strength. In order to find an accurate solution for both the FS and RS 
dynamics, we make use of Ull with the mechanical model. As explained above, this allows 
for the blast wave with a long-lived RS to satisfy the energy conservation, by introducing a 
pressure-gradient across the blast wave region. 

We perform a sophisticated calculation of afterglow emission, invoking a Lagrangian 
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description of the blast wave. In the widely-used analytical afterglow model (e.g. Sari et al. 
1998), it is assumed that the entire shocked material forms a single zone with same energy 
density and magnetic field. The electron energy distribution is solved only in the energy 
space, with no consideration of spatial distribution within the shocked region. Beloborodov 
(2005) (hereafter B05) described a more sophisticated Lagrangian method, in which the 
postshock region is resolved into subshells using a Lagrangian mass coordinate. B05 studies 
an evolution of the magnetic field and power-law spectrum of electrons for each subshell as 
the blast wave propagates. However, the postshock material in B05 is not resolved in radius 
and all the subshells are located at the same radius. Also, the pressure and energy density 
in B05 are assumed to be constant throughout the postshock material. Improving on B05, 
here we have a spatial resolution into the blast region, allowing for our Lagrangian shells to 
have their own radius. We further introduce a pressure profile that smoothly varies over the 
blast. This is because the pressure at the FS differs from the pressure at the RS, as discussed 
above. As the blast wave propagates, we keep track of an evolution of the pressure, energy 
density, and adiabatic index of every shell on the blast. We also keep track of an evolution 
of the magnetic field and power-law spectra of electrons of all shells. 

Finally, in order to calculate synchrotron radiation from a spherical shell on the blast, 
we analytically find an observed spectral flux for a distant observer, taking into account the 
effects of the shell's radial velocity and spherical curvature. We integrate this flux over the 
entire blast to find the sum of emissions from all the shells between the FS and RS. 

In Section [2J we briefly summarize how we find the dynamics of a blast wave with a 
long-lived RS. In Section [3j we describe in detail our method of calculating afterglow light 
curves. Numerical examples are presented in Section HJ which exhibit various features on 
the blast wave dynamics and afterglow light curves. Our results are summarized in Section [5] 
(Discussion) and Section E] (Conclusion). 

2. Dynamics of a blast wave with a long-lived RS 

A schematic diagram of a spherical blast wave is shown in Figure [TJ As mentioned 
in Section [TJ the RS wave is expected to be long-lived with a stratification on the ejecta 
shells. As the blast wave with a long-lived RS is not in the self-similar stage of BM76, its 
deceleration deviates from the solution of BM76. 

In order to find such dynamics of a blast wave with a long-lived RS, we make use of 
the semi-analytic formulation presented in Ull. The formulation has three input functions 
Pi(r), L e j(r), and r ej -(r), which define the initial setup of the blast wave (see Figure [Q. The 
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ambient medium density pi(r) is allowed to take an arbitrary radial profile, where r is the 
radius measured from the central engine. The ejecta is completely specified by two other 
functions, i.e., its kinetic luminosity L e ^{r) and Lorentz factor r e j(r). Here r indicates an 
ejection time of the ejecta shells. 

When two functions L^(t) and r e j(r) are known, a continuity equation, which governs 
a radial spread-out and spherical expansion of a stratified ejecta, can be solved. This then 
yields an analytic solution for the ejecta density p e j of any r-shell at radius r (Ull, Section 
3.1), 



Pcj(r,r)- L " 1 



l- r - 



1: 



c (r 2 - 1)3/2 

Here c is the speed of light, and v e j(r) = c (1 — l/r 2 ) 1 / 2 is the velocity of r-shell. Equation 
(CD) is exact for a non-increasing profile of r e j(r); ^(t) = dT^/dr < 0. 

We then self-consistently find the path of the RS wave through this evolved ejecta (Ull, 
Section 3.3). When the RS wave is located at radius ry(t) at time t, we numerically determine 
which r-shell gets shocked by the RS, and name it as r r (t)-shell; the subscript r in r r and r r 
refers to the RS. In other words, the RS wave sweeps up the r r (t)-shell at time t, which has 
the Lorentz factor r c j(r r ) = r c j(RS) and the density p e i{T r ,r r ) = p e j(RS). 

The formulation finds a dynamical evolution of the blast wave using the "mechanical 
model" (Beloborodov & Uhm 2006; Ull, Section 4). The mechanical model was developed 
for a relativistic blast wave by applying the conservation laws of energy-momentum tensor 
and mass flux on the blast between the FS and RS. Specifically, we numerically solve a set 
of coupled differential equations (78)-(80) and (92) of Ull, which makes use of the FS and 
RS jump conditions (Ull, Section 3.2). 

The FS and RS jump conditions in Ull are derived adopting a realistic equation of 
state (EOS) with a variable index k, 

p = K(e-pc 2 ), k=1m+IJ > (2) 

where p, e, and p are pressure, energy density, and mass density of the shocked gas, respec- 
tively, and 7 is the mean Lorentz factor of gas particles. The quantity k smoothly varies 
between 2/3 (for a non- relativistic gas) and 1/3 (for an ultra-relativistic gas) as the gas 
temperature varies. This EOS differs from the exact EOS of a relativistic ideal gas (e.g., 
Synge 1957) by less than 5 % (Ull, Section 2.2). Since e = 7pc 2 , the EOS in Equation ([2]) 
is the same as 
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which was previously introduced by Mathews (1971) considering a relativistic "monoener- 
getic" gas, and later adopted by Meliani et al. (2004) and Mignone & McKinney (2007) in 
their numerical simulations. 



3. Light curves from a blast wave with a long-lived RS 

Numerically solving for the dynamics of a blast wave with a long-lived RS, we discretize 
the ambient medium and ejecta into spherical mass shells 5m. At every calculation step, a 
pair of shells is impulsively heated by shock fronts; one shell by the FS and the other by the 
RS. We follow these shells subsequently and use them as our Lagrangian shells for the blast. 
Thus, the blast is viewed as being made of many different hot shells that pile up from the 
FS and RS, as depicted in Figure [1] with dotted (red) curves. 

Assuming an adiabatic blast wave, we keep track of an adiabatic evolution of these 
shells; an evolution of the thermodynamic quantities of shocked gas (pressure, energy density, 
adiabatic index, etc) is followed individually for every shell. This in turn yields an evolution 
of the magnetic field for the shell (Section I3.ip . We also keep track of an evolution of the 
power-law spectrum of electrons; a radiative and adiabatic cooling of the spectrum is followed 
for every shell (Section 13. 2p . 

Zooming in on the blast shown in Figure HJ we focus on a single spherical shell of radius 
r expanding with a Lorentz factor T (see Figure [2]). While taking into account the effects 
of the shell's Doppler boosting and spherical curvature on the synchrotron photons emitted 
along an observer's line of sight, we analytically find an observed spectral flux in terms of an 
observed frequency z/ obs and observer time t ^ s (Section l3.3p . As the blast is made of many 
Lagrangian shells, this spectral flux is summed over the blast; a sum of emissions from all the 
shells in the shocked ambient medium (or the shocked ejecta) is denoted by "FS emission" 
(or "RS emission"). 



3.1. Adiabatic evolution of the blast 

Here we follow the adiabatic evolution of the shocked gas on a shell, in order to find the 
evolution of the magnetic field for the shell. An adiabatic process of a relativistic gas whose 
EOS is specified by Equation ([2]) or fl3]) is described by (Ull, Section 3.4.1) 



— — = C = const. where VmW) = , tu • 

Pm{K) {l^-V 
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The function p m (/«) is monotonically decreasing in its valid range, ^ < k < 



3" 



Let us consider a shell 5m* either in the ambient medium or in the unshocked ejecta; 
an index i is added to specify this shell. When the shell 5m % is shocked by the FS (or 
the RS), the jump conditions of the FS (or the RS) determine its initial thermodynamic 
quantities: pressure p l , energy density e l , mass density p l , and quantity k 1 (Ull, Section 
3.2). Substituting the initial values p l and n l into Equation (j4]), we determine the constant 
C l of the shell Sm\ 

f -r = -A^. (5) 



Thus, if we know p l or k % of the shell 5m l at later times, we can subsequently follow an 
adiabatic evolution of the shell. 

Solving the mechanical model, at every calculation step, say at time t, we know Tf(t), 
r r (t), Pf(t), p r (t), and P(t), which are radii of the FS and RS, pressures at the FS and RS, 
and integrated pressure over the blast, respectively. Thus, an instantaneous pressure profile 
for the blast (r r < r < rf) may be approximated by a quadratic function p(r), which (1) 
matches two boundary values (i.e., p(rf(t)) = Pf(t) and p(r r (t)) = p r (t)) and (2) satisfies 
the integrated pressure P(t) = J^ f p(r) dr (Ull, Section 5). 

As we also know the radius r l of the shell Sm l at time t, we now have the pressure p l of 
the shell 5m l at time t: p l {t) = p(r l (t)). Equation ([5]) then allows us to numerically find the 
quantity k 1 of the shell 5m l at time t. All other thermodynamic quantities of the shell can 
be found accordingly. For instance, Equation (T5]) gives the thermal energy density e\ h of the 
shell as 

At) 

The electrons in the shell 5m l emit synchrotron radiation in a magnetic field B % . The field 
is unknown and parameterized by a micophysics parameter e#, which is the ratio of field 
energy density to thermal energy density: e# = (-^B l2 ) / e\ h . Thus, the magnetic field B l of 
the shell 5m l at time t is given by 

1/2 



l b (t)=e\t)-p\t)c> = ^. (6) 



B\t)= [87re B el h (t)] 



1/2 



, At) 
B W)\ 



(7) 



Hereafter, we will omit the index i to simplify our notation. 
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3.2. Power-law spectrum of electrons 

We assume that a non-thermal electron spectrum is created in a fresh shell 5m at a shock 
front (FS or RS); i.e., the electrons are accelerated into a power-law distribution above a 
minimum Lorentz factor 7m , 

dN 

/( 7e ) = _= J K'( 7e -l)-P for 7e > 7m , (8) 

where K is a constant, p is the slope of the spectrum, and 7e is the Lorentz factor of the 
accelerated electrons in the fluid frame. The total number of electrons within the spectrum 
is given by 

SN= f(l e )d le = _( 7m -l)i-P. (9) 

hm p - 1 

The thermal energy of all electrons in the spectrum is found as 

f°° Km r 2 

SE = / [( 7e - l)m e c 2 ] /( 7e ) d le = ( 7m - l) 2 ^, p > 2, (10) 

where m e is the electron mass. Thus we find the mean thermal energy per electron, 

5E p-1 



5N p-2 



( 7m - 1) m e c . (11) 



For neutral plasma without pair-loading, if (1) all electrons passing through the shock become 
non-thermaU) and (2) a fraction e e of the shock energy goes to the electrons^], then the mean 
thermal energy per electron is alternatively given by 

8E ._ . n 

= e e {lp - 1) m p c . (12) 

Here m p is the proton mass and 7p is the mean Lorentz factor of the protons in the postshock 
medium. Equations (iTTj) and (fT2|) are combined to yield the lowest Lorentz factor 7m as 

p — 2 m v , , . 

7m = l + ^— r— e e (7 P -l ■ 13 
p — 1 m e 

For a fresh shell created at a shock front, 7p equals the shock strength, i.e., the relative 
Lorentz factor of the preshock to the postshock medium (BM76; Ull, Equation (14)). Thus, 



1 This is a usual assumption in the afterglow literature, although it is not necessarily true. See Genet et 
al. (2007) for examples of afterglow light curves in the long-lived RS scenario, obtained by assuming that 
only a small fraction £ of the electrons become non-thermal. 

2 The dominant fraction of shock energy goes to the protons, which dominate the pressure of the blast 
and evolve adiabatically. 
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for the fresh shell at the FS, 7 P equals T, the Lorentz factor of the blast, as the ambient 
medium is at rest in the lab. frame. For the fresh shell at the RS, 7p equals 743 given by 

74 3 = rr CJ -[(r 2 -i)(rJ-i)] 1/2 , (14) 

where r e j is to be evaluated for the r r -shell; r e j = r e j(RS). 

We use Equation (J5J) to find the constant K and substitute it into the electron spectrum, 



f{le) = 5N P 



7n 



(15) 



7m 

This spectrum created at the shell 5m evolves as the blast wave propagates. We track the 
adiabatic and radiative cooling of the electron spectrum as follows. 

1. Adiabatic cooling of 7 m : an adiabatic cooling of relativistic electrons is described as 
7 e oc p 1 / 4 . Thus the minimum Lorentz factor "f m also evolves as 7 m oc p 1//4 . Here p is the 
pressure of the shell, not the slope of the electron spectrum. 

2. Radiative and adiabatic cooling at high 7 e : the electrons at high j e in the spectrum 
experience a radiative and adiabatic cooling, which is described by the first and the second 
term below, respectively, 

7e = -^— 5 2 (l + >07 e 2 + ^7e. (16) 

D7T m e c 4 p 

Here <jt is the Thomson cross section, and the dot indicates a derivative with respect to t', 
the time measured in the co-moving fluid frame. The Compton parameter Y describes a 
relative contribution of inverse Compton scattering to the cooling rate of electrons. 

Equation ([TBI) defines a cutoff Lorentz factor 7 C at the high end of the spectrum. Dividing 
Equation ([16]) by — 7J, we get 

d f 1 \ 1 a T 2 M i ^ 1 f l \ dlnp 



B 2 (l + Y)-- - —f. (17) 



dt' \7c/ 67r m e c 4 \7 c y dt 

The parameter Y can be evaluated as (Sari & Esin 2001 



where 



(l + Y)Y = ^ V or Y = -l + J\ + ^ V , (18) 

v = i (im/lcf 2 for 7 m <7 c , 

I 1 for 7 C < 7 m . 

Note that p in Equation (TTT]) is the pressure of the shell, and p in Equation (Tl9~]) is the slope 
of the electron spectrum. When a fresh shell 5m is created at a shock front (FS or RS), 
7 C = +00 is adopted as its initial value. We then solve Equation f JT7|) numerically and find 
a subsequent evolution of 1/7 C (and hence, 7 C ) for the shell 5m. 
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3.3. Curvature effect and light curves 

Consider a spherical shell 5m on the blast, which has radius r at time t expanding with a 
Lorentz factor T. An observer is located in the positive ^-direction at a large (cosmological) 
distance. See Figure [2] for a schematic diagram. An observer time t ^ s is set equal to zero 
when the observer detects the very first photon that was emitted at the explosion center at 
time t — 0. Also, consider a thin ring on the shell 5m, which is specified by a polar angle 9 
with respect to the z-axis. Then photons emitted from this ring at time t will be detected 
by the observer at the observer time, 



t 



obs 



r 

t ji 

c . 



(1+z), (20) 



where fi = cos 9 and z is the cosmological redshift of the burst. The factor (1 + z) is 
introduced due to the expansion of the universe. 

The spherical shell 5m contains a total of 5N electrons uniformly distributed over the 
shell, retaining a non-thermal spectrum of j m , j c , and slope p (see Section I3T21) . Since the 
number of electrons on the thin ring (between 9 and 9 + 59) is given as (\5fj,\/2) 5N, the 
electron spectrum of the ring becomes 



~ \5fi\5N p-l 

/(7e) = —r~ V" 



In. 



(21) 



where we assumed 7 e 3> 1 and 7 m 3> 1. Here the tilde on /(7 e ) indicates the ring. Note 
that /(7e) is subject to the thickness of the ring, 59 or 5/i. 

We now focus on a single electron in the ring. On average, the electron is assumed to 
emit synchrotron photons isotropicalhjf] at a single characteristic frequency vauid in the fluid 
frame, 

^uid = (0.15)^7e- (22) 

m e c 

Here q e is the electric charge of the electron, and B is the magnetic field of the shell. The 
synchrotron luminosity of the electron is given by 

L e = ia T cB 2 i (23) 

The photon frequency z/j a b in the lab. frame is different from z/fl ui d due to radial expansion 
of the shell. Since the radial bulk motion of the electron has the angle 9 with the observer's 



3 See Beloborodov et al. (2011) for discussion of anisotropic emission in the fluid frame. 
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line of sight (see Figure [2]), the frequency v\ a \> in the direction of the observer is simply 

r(l — p cos 9) 

The angular distribution of the energy of photons emitted by the electron in the direction 
of the observer is written in the lab. frame as 

dE 1 \ L e 1 (25) 



dndtj lah 4vr T 4 (l -/3cosfl) 3 
Here the subscript in E\ refers to the single electron. 

From Equations (122 p and (1241) , we note that the Lorentz factors 7 m and 7 C of the electron 
spectrum correspond to the frequencies z/^f 3 and z/£ , respectively, 

^-^ q £w^-y (26) 

The electron with 7 e = 7^ lab has the frequency i^ a b in the direction of the observer, 

^ = (0.15)^ — ^ (28) 
m e c 1 (1 — p/jL) 

The spherical shell emits photons continuously when it expands, but we think of a 
series of "snapshots" of the shell. We consider two consecutive snapshots (i.e., two calculation 
steps) separated by a time interval St, and then assume that the shell accumulates its emission 
between two snapshots and emits all the accumulated energy instantaneously like a "flash" 
when it arrives at the second snapshot. In other words, the emission from the shell is viewed 
as a series of flashes. 

When the shell flashes, the accumulated energy by the single electron during 5t is 
emitted. The energy emitted into a solid angle 5Q in the direction of the observer is given 
in the lab. frame as 

The next step is then to calculate the emission from the entire thin ring; let 5E lah be the 
emission from the entire ring into the solid angle 5Q in the direction of the observer. With 
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the definition of spectral energy 5El ah = d(5E l&h )/dv } we consider (e.g., B05) 



lab \ 
v\ah> 



d(5E 



lab\ 



d(ln v) 



d(5E 



lab \ 



"="lab 

lab"\ 



2rf(ln 7e 



7e=7i' 



7"l a b 



d(5E lab ) 



dje 



(30) 



7e=7», s 



Here the second equality uses the relation 5(lni^ a b) = 2 5(ln7 nab ), which is verified from 



Equation (1281) ; for the instantaneous flash of the shell, a variation in v^b results only from 
5 7l /i ab since the ring has a fixed B, T, and /x instantaneously. Together with Equations (j2ip 
and d29]), Equation (00]) yields 



^ab(^) 



7"lab 



1 P 



247T 47T 



otc -E> 2 7 , 2 



"lab 



Te=7 nab 

5Q5t \Sfj,\SN 
r 4 (i - M 3 



. 7m 



1 1-p 



(31) 



We remark that Equation ( l3Tj) is valid only for z/^ b < z/i a b < since the electron spectrum 
in Equation ( l2~Tj) is valid only for 7 m < 7^ ab < 7c . Dividing Equation (l3*T]) by z/i ab in Equation 
( 1281) . we find the spectral energy emitted from the entire thin ring into the solid angle 5Q in 
the direction of the observer, 



5E 



lab 

"lab 



5 p- 1 cr T m e c 2 B 5Q5t \5fi\5N 
q e r3(l-/3/x) 2 



18vr 4tt 



T^la 

7m 



1-p 



(32) 



This energy is emitted instantaneously during the flash, but the thickness of the ring (between 
9 and 9 + 59) introduces a time interval 5t\ ah along the observer's line of sight, 



^iab = - [cos 9- cos(6> + 89)} = -\6fj,\. 
c c 



(33) 



Hence, Equations (I3"2"j) and ( 13"3"|) yield the spectral luminosity of the entire ring, 



5L 



lab 

"lab 



5E)f° 
"lab 

<^lab 

5 p - 1 <7 T m e c 3 5 <5fi 5t <5iV 



187T 47T 



g e r r 3 (l-/3/i)^ 



T^a 

7m 



1-p 



(34) 



which shines into the solid angle 5Q in the direction of the observer. Note that we do not 
include a tilde for 5L^ since the thickness of the ring (namely, 5 fx) cancels out here. 
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The burst is at a cosmological distance from the observer. For simplicity, we consider a 
flat matter-dominated universe. The luminosity distance of a burst at redshift z is 



Dr = — 



2c 
Ho 



1 + z)- y/(l + Z 
(2.6 x 10 28 cm) [(1 + z)- y/(l + z) 



(35) 



where Hq ~ 70 km s 1 Mpc 1 is the Hubble constant. The photons are redshifted while 
traveling the cosmological distance, and the observed frequency z/ obs is obtained by 

^lab 



^obs 



l + z 



(36) 



The definition of the luminosity distance gives an observed spectral flux at z/ bs, 



SF; 



-obs 



5(1 + z) p-1 o T m e c 3 B 5N St 



^lab 



1 (1-J0/2 



lab 



18tt AttDI q e r r 3 (l-/3/i) 2 
where Equations f )26|) and fl28|) have been used to replace 7^ ab /7m by (^lab/^m 13 ) 1 ^ 2 - 



(37) 



Finding // = (c/r) [£ — t Q bs/ (1 + z )] from Equation fj20|) and substituting it into Equation 
( I37p . we arrive at an analytical expression for SF° h ^ in terms of t Q b s and z/ Q b s , 



SF^(t ohs ,u ohs ) = SF l 



^obs 



X 



^obs 
,,obs 
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where 



5F max 



,obs 



5(1 + z) p — 1 <Jt m e c 3 B 



SN St 



18tt AttD'I 
0.15 g e 5 



' (71 



r 3 [i-f (*-!£)] 



2 ' 



1 + 2 m e c r [l - ^ (t - f^)] 



(38) 

(39) 
(40) 



When the spherical shell of radius r flashes at time t, the observer at cosmological distance 
receives its emission for a period of the observer time t Q b s , 



r 

£- - < 

c ~ l + z 



f,,,,s <t+~. 

c 



(41) 



In other words, Equation (j38j) is to be evaluated only in this period of t Q b s - For any t Q b s 
in this period, the observed spectral flux SF°^ a at the observed frequency z/ obs is expressed 
analytically. Thus we do not need to consider the ring any more. 
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Equation fl3T|) is valid for z/^ b < < z/^ ab , and therefore Equation fl38l) is valid only 

for z/° bs < v ohs < z/° bs , where 



,„,,s _ 0-15 g e £ t| 

1 + 2 m e c r [l - & (* - |^)] 



^ = T— : TTT ^ rrvT - ( 42 ) 



In general, the observed spectral flux at z/ Q b s can be obtained as 

' (^obs/C bS ) 1/3 for Z/ obs < z/f < ^° bs 

(^obs/C bs ) (1 " P)/2 for C bs <^obs<^ 



^obs ^obs 



C 1 

(^obs/^ c obs ) 1/3 for z/ obs < v° hs < z^ bs , 

for z/° bs < z/ obs . 



(43) 



Equation ( 1431) needs to be summed over the shocked region as the blast is made of many 
different shells. 

Here we recover the index i to name our Lagrangian shells {5m*}. Each 5m l is impul- 
sively heated at some point by the FS or the RS and becomes a shocked shell on the blast. 
We denote these shocked shells by {Sml hd }. As mentioned above, the shells {Sml hd } on the 
blast have their own individual radius r l , number of electrons SN l , magnetic field B % , and 
the Lorentz factors Y m and 7*. 

Let us now use another index j to specify the time P of each calculation step (or flash). 
Solving for the blast wave dynamics at a calculation step with time P , we find the Lorentz 
factor T(P) of the blast and the radii {^(P)} of the shells {Sm l shd }. We also evaluate the 
shells' emission properties {B l (P), Y m (P), 7* (P)} at time P (Sections 13. II and I3.2j) . 

Then, Equation fT43|) gives the spectral flux <5i™ b ° bs of each shell 5m* hd at time P . We 
sum it over all shocked shells \Smt hr] } to find a flux 5F?. obs = Y] IXm i 1 5FJ? obs of the entire 
blast at time P . Lastly, this needs to be summed over all flashes separated by time intervals 
{SP} to give a total spectral flux, 

^(tobs^obs) = E E ^ b f S (Us^obs). (44) 
{SP} {5ml hi } 

Note that two summations are not commutative since the FS and RS waves create new 
shocked shells as time goes. We find the "FS emission" (or the "RS emission") by taking 
the summation Y^{Srri } over an snens only in the shocked ambient medium (or only in 
the shocked ejecta). For a fixed observed frequency i/ b s , Equation (H4|) gives light curves 
F° o bs s (t D bs) at i/ bs- For a fixed observer time t b s , Equation (jUJ) yields flux spectra F° o bs s (z/ b s ) 

&t ^obs- 
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4. Numerical examples 

Now we present a total of 20 different numerical examples, which are named as follows; 
la, lb, lc, 2a, 2b, 2c, 3a, 3b, 3c, 4a, 4b, 4c, 4d, 5a, 5b, 6a, 6b, 6c, 6d, and 6e. For all 20 
examples, we keep the followings to be the same: (1) A constant density pi(r)/m p = 1 cm -3 
is assumed for the ambient mediunfl (2) The ejecta has a constant kinetic luminosity 
L C j(r) = L = 10 53 erg/s for a duration of r b = 10 s, so that the total isotropic energy of the 
burst is to be E b = L r b = 10 54 ergs. (3) The burst is assumed to be located at a redshift 
z = 1. (4) The emission parameters are e e = 10 _1 , €b = 10~ 2 , and p = 2.3 for the RS light 
curves, and e e = 10 -2 , e# = 10 -4 , and p = 2.3 for the FS light curves. 

Note that we have adopted different micophysics parameters for the FS and RS. This 
is because the FS and RS shocked regions originate from different sources and the strengths 
of the two shocks can be significantly different. Indeed, afterglow modeling suggested that 
the RS can be more magnetized, and e e of the two shocks can also be different (e.g. Fan 
et al. 2002; Zhang et al. 2003; Kumar & Panaitescu 2003). Since the GRB central engine 
is likely magnetized, a natural consequence would be to invoke a larger in the RS. The 
bright optical flash observed in GRB 990123 (Akerlof et al. 1999) is likely related to such a 
case (Zhang et al. 2003). Our emission parameters above are chosen such that the FS and 
RS fluxes are comparable to each other. Of course, either the FS or the RS emission could 
be further enhanced or suppressed by varying their parameters e# and/or e e . 

The only difference among examples then goes on the profile of the ejecta Lorentz factors 
r c j(r) as a function of ejection time r; i.e., the examples have a different ejecta stratification 
within the same duration r b = 10 s. Since we are mainly interested in the afterglow light 
curves, we ignore the initial variation of r e j and take a simple uniform profile of high Lorentz 
factors early on. Thus, for all 20 examples, we assume a constant Lorentz factor r e j = 500 
for the initial 3 s. From 3 to 10 s, the examples have a decreasing profile of r e jjf| exhibiting 
various types of the ejecta stratifications. Note that only a comparable amount, 70 % of 
the burst energy, has been distributed over the shells with lower Lorentz factors in order to 
maintain a long-lived RS. Thus, the deceleration of the blast wave would deviate only mildly 
from the solution of BM76. 



4 Although our method described in this paper allows us to study other types of environment, such as a 
stellar wind with pi(r) oc r~ 2 , we will keep the same density p\{r)/m v — 1 cm" 3 in all 20 examples. This is 
to focus on the study of the ejecta stratifications. 

5 It is natural to assume a decreasing profile of T e j since the internal shocks during the prompt emission 
phase tend to smoothen the r c j distribution and lead to a decreasing r c j profile (otherwise, additional internal 
shocks would occur). 
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The examples with the same number in their names share a similar shape of the ejecta 
stratifications, and we categorize 20 examples into 7 different groups; (la/lb/lc), (2a/2b/2c), 
(3a/3b/3c), (4a/4b/4c/4d), (la/5a/5b), (6a/6b/4d), and (6c/6d/6e). This is to provide an 
efficient comparison among examples. Note that we use la and 4d twice in the comparisons. 

For each group of examples, we present 3 figures: (1) The 1st figure shows the ejecta 
stratifications (e.g., see Figure |3]). The ejecta Lorentz factors r e j(r) are shown as a function 
of ejection time r, in different line (color) types. (2) The 2nd figure shows the blast wave 
dynamics of the examples (e.g., see Figure H]). In the panel (a), the curves denoted by T show 
the Lorentz factor of the blast wave as a function of the radius r r of the RS, and the curves 
by r e j(RS) show the Lorentz factor of the ejecta shell that gets shocked by the RS when the 
RS is located at the radius r r . In the panel (b), the curves denoted by pf show the pressure 
at the FS as a function of the radius r r , and the curves by p r show the pressure at the RS. 
Since the ambient medium is assumed to have a constant density pi, the pj curves resemble 
the T curves; pf oc T 2 pi. The panel (c) shows the relative Lorentz factor 743 across the RS 
wave, as is given by Equation ( ff4"]) . The panel (d) shows the density n e j(RS) = p C j(RS)/m p 
of the ejecta shell, which enters the RS wave at radius r r . We omit the panels (c) and (d) for 
the last two groups of examples. (3) The 3rd figure shows the afterglow light curves (e.g., see 
Figure |5]). In the upper panel (a), we show the FS emissions in X-ray (1 keV) and R band 
as a function of the observer time £ t, s . In the lower panel (b), we show the RS emissions in 
X-ray (1 keV) and R band. 

Each group has a comparison point to help readers understand the effects of ejecta 
stratification on the FS and RS dynamics and further on the FS and RS afterglow light 
curves. 

4.1. Group (la/lb/lc) 

Figure [3] shows the ejecta stratifications of examples la, lb, and lc. For a duration from 
3 to 10 s, the ejecta Lorentz factors r e j(r) decrease exponentially from 500 to 5, y/5 x 50, 
and 50, respectively. Note that an exponential decrease implies jj:(lnr c j) = Fgj/r e j = const. 

The blast wave dynamics of these three examples are shown in Figure HI The r e j(RS), 
p r , 743, and pej(RS) curves vanish at an earlier time or radius for higher ending values of r ej 
(i.e., examples lb and lc), as the RS waves cross the end of ejecta. The T curve of example la 

3/9 

shows that its blast wave decelerates slightly slower than T cx r r , the self-similar solution 
of BM76 (denoted by a dot-dashed line). The examples lb and lc have higher RS pressure p r 
than the example la, and therefore their T curves deviate from BM76 even stronger than the 
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case of la while their RS waves exist. However, once their RS waves vanish, their T curves 
start to follow T oc r^ 3 / 2 as they should. As a result, all three T curves completely agree with 
one another after all three RS waves disappear. This is expected because the same amount 
of burst energy Ef, has been injected into the blast waves which swept up the same amount 
of ambient medium out to a certain radius. In other words, three T curves are shaped by 
three different plans or time schedules of "spending" the same energy budget E\,. Once the 
budget is used up, the outcome or the Lorentz factor of the blast waves should be the same. 
Also notice that three pf curves exhibit the same behaviors as V curves, accordingly, since 
Pf oc T 2 pi and pi = const. 

Three r e j(RS) curves stay close to the T curves since the ejecta shells catch up with 
the blast waves only when r e j ~ T; the resulting relative Lorentz factors 743 are shown in 
panel (c). A constant kinetic luminosity L e ^(r) = L is assumed here, and thus Equation 
(TT]) yields p e j(r, r) oc r -2 when Tgj(r) = 0. Hence, p e j(RS) = p E j(T r ,r r ) oc r~ 2 while the RS 
waves sweep up the initial 3 s of ejecta shells with r e j = 500 (see panel(d)). When the RS 
waves arrive at r = 3 s, they encounter a discontinuity in the value of T'^(t), which results 
in a sudden drop in p e j(r, r) across r = 3 s; see Equation ([T]). Note that both p e j(RS) and 
p r curves exhibit a sudden drop correspondingly. 

When Tgj(r) < 0, Equation fll}) simplifies 

p ej (r,r) oc M^0 oj ( r ) if g ci {r) < (45) 
r ci ej 

where we have defined 

-1 

(46) 

For exponentially decreasing parts of r e j (r > 3 s), the examples la, lb, and lc have a 
constant value for (7 C j( r ), which equals 1.52, 2.03, and 3.04, respectively. Since -^r ~ t Q b s , 
Equation fT45]) is applicable here, and yields p e j(T, r) oc r -3 . Hence, p e j(RS) = p e j(T r , rv) oc r~ 3 
while the RS waves sweep up these exponential parts of ejecta shells (see panel (d)jj. Then, 
the p r curves roughly follow p r oc r~ 3 , as the 743 curves mildly increase at around the value 
1.1. Note that the pf curves also roughly follow p/ oc r^ 3 , since the T curves deviate only 
mildly from T oc r^ 3 ^ 2 . 

Therefore, it is not surprising that both the FS and RS light curves exhibit power-law 
declines with a very similar value of temporal indices (see Figure [5]). Note that our choice 



6 It is evident here that a radial spread-out of a stratified ejecta induces a lower density ejecta-flow than 
a non-stratified ejecta with r c j = const, due to its gradient r^(r). 



T 
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of emission parameters e e and e# produces the RS afterglow emissions at a comparable flux 
level to the FS emissions. An abrupt decline of the RS light curves after a rising phase is 
due to a sudden drop in the p r curves (mentioned above); this will become clear below with 
the group of examples la, 5a, and 5b. 

Consequences of spending the energy budget in three different ways can be evidently 
seen in both the FS and RS light curves. In particular, the RS waves of examples lb and lc 
vanish at an earlier time or radius, and thus their RS light curves exhibit a temporal break 
and steepen afterwards. Even after the RS waves disappear, the RS light curves still get 
contributions from previously shocked Lagrangian shells whose z/° bs (Equation ( H2|) ) is still 
higher than z/ b s ; namely, the RS light curves are produced by residual emissions plus high 
latitude emissions. In the case of examples lb and lc, the RS light curves are essentially 
governed by high latitude emissions after the temporal break. Notice that this RS break 
looks like a "jet break", which has been interpreted to be caused by a collimation of ejected 
outflow. 



4.2. Group (2a/2b/2c) 

The ejecta stratifications of examples 2a, 2b, and 2c are shown in Figure O After a 
steep decrease at r = 3 s, the r e j profiles transition to an exponential decrease, displaying 
three different concave shapes. 

A large gradient r^(r) at r = 3 + e (e < 1) results in a very small value for g C ](j). 
The g C )(j) then gradually increases as it goes through the concave shapes, and becomes a 
constant value of 3.26 when it continues on the exponential parts. Thus, Equation (I45p is 
applicable here. Hence, p e j(RS) = p${T r ,r r ) oc r~ 3 g e $(r r ). Note that, when g e j(T r ) rises, it 
competes with a decaying term r~ 3 to determine the density profile p e j(RS) of ejecta shells 
entering the RS waves. As shown in Figured the p C j(RS) curves exhibit a rising phase and 
then eventually follow p e j(RS) oc r~ 3 when the RS waves sweep up the exponential parts 
of ejecta shells. A local maximum of the p e j(RS) curves is located at different radii r r as a 
result of the competition between g e ]{j r ) and r~ 3 . 

The p r curves resemble the p e j (RS) curves; a significant drop and then a gradual recovery, 
followed by a power-law decline, p r oc r~ 3 roughly (see panel (b)). While the RS pressures 
p r are weakened temporarily, the blast waves satisfy p r <C P/, and therefore the T curves 
follow r oc TV . When the p r curves start to recover from the drop, the V curves start to 
deviate from the self-similar solution. Note that three T curves show a very minor difference 
because a relatively large difference in p r curves during the recovery phase does not have 
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much meaning to the T curves since p r <C Pf. Three pj curves also show a negligible difference 
accordingly. 

Hence, the FS light curves do not exhibit any noticeable difference (see FigureED- On the 
other hand, three different p r curves, of course, show up in the RS light curves. During the 
temporary weakening of the RS waves, their X-ray light curves are dominated by high latitude 
emissions since the synchrotron frequencies v c corresponding to the cutoff Lorentz factors 7 C 
of previously shocked Lagrangian shells are already below 1 keV. However, the frequencies 
v c stay above the R band for some time for the given values of emission parameters; recall 
that the Lorentz factor j c of each Lagrangian shell is determined by Equation ffTT|) . which 
involves e e , €b, p, etc. Therefore, even when the RS waves have been weakened significantly, 
the residual emissions from these previously shocked shells still contribute to the R band, 
and produce the R band light curves with a shallower decline (temporal index a ~ 1) than 
the X-ray light curves of high latitude emissions (a ~ 3)0. 

During the recovery phase of the RS waves, their increasing p r curves are combined 
with decreasing T curves, and as a consequence, the X-ray light curves exhibit a flattening 
phase ("plateau"). Note that the detailed shape of this plateau phase is inherited from the 
p r curves that are shaped by the concave parts of r e j profiles (Figure [6]). This recovery of 
the RS waves is not fully shown up in the R band light curves since it is buried under the 
residual emissions (mentioned above). For a smaller value of e# in the RS, the electrons 
in the shocked shells would cool more slowly, and their resulting residual emissions would 
continue for a longer period of time in the R band, thus enhancing the chromaticity between 
the X-ray and R band. 

As the RS waves continue on the exponential parts of r e j profiles, the p r curves roughly 
follow p r oc r~ 3 , and the RS light curves show the usual power-law decline with a temporal 
index (a ~ 1), very similar to that of the FS light curves. A temporal break at the end of 
the X-ray plateau phase marks the end of a rising phase of the p r curves, and therefore this 
"dynamically- caused" temporal break does not involve a spectral evolution across the break. 

4.3. Group (3a/3b/3c) 

The r e j profiles of examples 3a, 3b, and 3c are composed of a steep drop and a concave 
part, followed by a convex part (see Figure [H])- The example 3b has a higher ending value of 



7 Note that the high latitude emission from an end of the prompt internal shocks would also naturally 
produce a similar steep decline in X-ray (Hascoet et al. 2012). 
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The g e j(T r ) functions increase until the end of concave parts (i.e., inflection points), and 
then decrease afterwards through convex parts. Thus, the p C j(RS) curves rise, exhibit a 
local maximum, and then decay afterwards (see Figure [TOl . Since p e j(RS) oc r~ 3 g C j(r r ), the 
location of this local maximum does not exactly agree with, but roughly coincides with, the 
inflection point of lnr e j. Since the functions g e j(T r ) decrease after the inflection point, the 
Pej(RS) curves exhibit a more pronounced local maximum and decay faster afterwards than 
the case of examples 2a, 2b, and 2c. 

In particular, the example 3c shows a well-pronounced local maximum in the p e j(RS) 
and p r curves (Figure [TOl . which appears in the RS light curves as a small flaring activity at 
the end of the plateau phase (Figure [TT]) . The examples 3a and 3b also show a clear temporal 
break at the end of the X-ray plateau phase. As the RS wave of example 3b vanishes at an 
earlier time or radius, its RS light curves display an additional temporal break at t b s ~ 10 5 s. 



4.4. Group (4a/4b/4c/4d) 

The ejecta stratifications of examples 4a, 4b, 4c, and 4d have the same concave shape 
with various ending values that are equally spaced (Figure [T2"j) . 

Initially, the T curves follow T oc r r when p r C pj, and then start to deviate from 
it as the p r curves recover from a weakening (Figure [T3]) . Soon afterwards, the RS waves 
vanish, and their T curves again follow T oc rr 3 ^ 2 . All four T curves, of course, agree with 
one another after the RS waves disappear. 

The pj curves resemble the T curves, and the FS light curves show corresponding be- 
haviors (Figure [T4l) . During the recovery phase of the p r curves, their X-ray light curves 
produce a plateau phase. This plateau phase ends as the RS waves cross the end of ejecta. 
Therefore, the X-ray light curves exhibit a sudden steep decline (a ~ 3) beyond this plateau 
phase. Note that, the less shallow the plateau phase is, the longer it is. While the example 
4d has an almost flat plateau phase (a ~ 0), the example 4a has a mild plateau (a ~ 0.5), 
which is so long that its end can not be seen here until t Q b s = 10 6 s. 



4.5. Group (la/5a/5b) 



New examples 5a and 5b are compared to a previous example la here (Figure [TBI) . The 
examples 5a and 5b show no steep drop at r = 3 s. Instead, their r e j profiles are made of a 
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single convex distribution. 

The example 5b has no discontinuity in r^(r) across r = 3 s, and therefore, its p e j(RS) 
curve shows no sudden drop at r = 3 s (Figure [TBI) . When F' e -(r) ~ beyond r = 3 s, 
Equation flTJ roughly yields p e j(r, r) oc r -2 , and its simplified version, Equation (145]) . is 
not valid; hence, p e j(RS) = p C j(r r ,r r ) oc r~ 2 . For the convex part of r e j profile, the g C j(r) 
function keeps decreasing as r increases, and as a result, Equation ( 145]) becomes applicable 
afterwards; p e j(RS) oc r~ 3 g e $(r r ). Thus, the p e j(RS) curve starts to decrease faster than that 
of example la, for which p e j(RS) oc r~ 3 . 

The p r curves resemble the p e j(RS) curves (panel (b)). Since the p r curve of example 
5a is initially higher than that of example la, its blast wave deceleration (i.e., F curve) 
is delayed more significantly than that of example la. The T curve of example 5b follows 
T oc ry 3 ^ 2 afterwards when p r <pj. 

The pf curves and the FS light curves show the same behaviors correspondingly (Fig- 
ured?]). The RS light curves also exhibit a similar behavior, but with a stronger variation in 
their temporal indices, since the p r curves have stronger variation than the pf curves. Note 
that the example 5b has no abrupt decline in its RS light curves after a rising phase since it 
has no sudden drop in its p r curve. This group of examples demonstrates that the FS and 
RS light curves can exhibit various temporal indices without varying the slope p of electron 
spectrum. In particular, this becomes significant at late times in the RS light curves. 

4.6. Group (6a/6b/4d) 

New examples 6a and 6b are compared to a previous example 4d (Figure ITS1) . Three 
examples have the same T e j profile up until r = 7 s. From r = 7 s, the example 6a has an 
exponential decrease while the example 6b continues with a convex profile. 

The p r curve of example 6a has another sudden drop at r = 7 s, since F' e -(r) is not 
continuous there (Figure [T9l). The p r curve then roughly follows p r oc r~ 3 as the RS wave 
sweeps up the exponential part of ejecta shells. On the other hand, the p r curve of example 
6b has no drop at r = 7 s, and then decays faster than p r oc r~ 3 , just like the example 5b. 

All three examples produce a flat plateau phase in the RS X-ray light curves while their 
RS waves sweep up the ejecta shells up until r = 7 s (Figure [201). The RS wave of example 
4d extends the plateau phase slightly as it sweeps the remaining ejecta shells (from 7 to 10 
s), and then vanishes. Thus, its X-ray light curve has a ~ 3 (high latitude emission) beyond 
the plateau phase. The RS wave of example 6b does not disappear, but becomes weaker and 
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weaker, producing its X-ray light curve with a ~ 2 beyond the plateau phase. Lastly, as the 
RS wave of example 6a continues on the exponential part, its X-ray light curve shows the 
usual a ~ 1 phase. A steeper decrease between the plateau and a ~ 1 phase is due to a 
sudden drop of p r curve at r = 7 s. 

4.7. Group (6c/6d/6e) 

Like the examples 6a and 6b, three examples 6c, 6d, and 6e share the same r e j profile 
with the example 4d up until r = 7 s (Figure l2~Tj) . From r = 7 s, three examples have a 
single concave shape with various ending values of r e j. 

The blast wave dynamics are shown in Figure [22j The RS light curves show a X-ray 
plateau and then a steep decay (a ~ 3), followed by three different types of re-brightenings 
(Figure [22}. 

5. Discussion 

The T e j profiles of all 20 examples are shown together in Figure [2H the line (color) type 
of each example is the same as previous sections. Each example plans on its schedule of 
spending the same energy budget Eb = 10 54 ergs, by shaping its r e j profile (i.e., its ejecta 
stratification). Depending on this schedule, its T curve goes through various deceleration 
path, but eventually arrives at the same spot in the T - r r plane. Indeed, the T curves of all 
20 examples arrive at the same spot as is shown in Figure [251 

The consistency shown in Figure [25] could serve as a "test" that a blast wave modelling 
with a long-lived RS needs to pass. While adopting a customary pressure balance p r = Pf 
across the blast wave, we find the blast wave dynamics again for examples la and 5b. The 
obtained T curves, denoted by la c and 5b c , respectively, are shown in Figure |26| together 
with previous T curves la and 5b (taken from Figure [TBI) . 

The T curves denoted by la c and 5b c clearly show that the pressure balance p r = Pf 
does not pass the test above. Moreover, a sudden drop or deceleration shown in the la c 
curve is not physical, since it is a mere consequence of enforcing p r = pj on the blast when 
the RS wave encounters a sudden drop in the p e j(RS) curve; recall that the example la has 
an abrupt drop in its ejecta density (Section 14. ip . The example 5b has no sudden change 
in its p e j(RS) curve (Section 14.51) . However, the 5b c curve also reveals a problem with the 
pressure balance p r = Pf, as it decreases faster than T oc r r 3 ^ 2 even though it is meant to 
describe an adiabatic blast wave for p\ = const. An existing long-lived RS would only delay 
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the deceleration of its blast wave, as shown by the la and 5b curves. Hence, this implies that 
the blast wave dynamics described by the 5b c curve does not satisfy the energy conservation. 
In fact, the la c curve does not conserve the energy, either; its arrival spot is very different 
from that of the la and 5b curves. 

Thus, enforcing an equality p r = pj or a constant ratio Pf/p r = const, on a blast wave 
with a long-lived RS would give rise to an incorrect dynamics. Furthermore, those methods 
would not be able to "capture" distinctive features of the RS dynamics (that are illustrated 
with 20 examples above), since the FS and the RS wave would force each other to follow a 
similar dynamics. 

With a mild energy injection into the blast wave, its deceleration (or the FS dynamics) 
deviates only mildly and smoothly from the solution of BM76, as is visible in Figure [25] 
In the meantime, its RS wave responds sensitively to the density flow p e j(RS), shaped by 
various ejecta stratifications, and produces fast and strong evolutions in the RS dynamics. 
These distinctive features then show up in the RS light curves, so that the RS light curves 
exhibit more diverse and vigorous behaviors than the FS light curves. This can be clearly 
seen in Figure [27] where we put together the light curves of all 20 examples. 

We remark that the diversity and non-trivial features shown in the RS light curves 
(Figure [27]) are produced as we simply move through the space of ejecta stratifications 
(Figure 121]) . All 20 examples have the same density pi(r)/m p = 1 cm~ 3 for the ambient 
medium, and the same duration = 10 s and isotropic energy = 10 54 ergs of the burst. 
The microphysics parameters are also kept to be the same in all examples. Thus, contrary 
to other proposed scenarios, this diversity does not require strong variations on the life time 
or the energetics of the central engine. 

In fact, many features shown in the RS light curves resemble what has been observed in 
the afterglow light curves. In particular, a weak bump at the end of the plateau phase, as in 
example 3c (Figure [TT]) . has been observed in the X-ray light curve of a few GRB afterglows 
such as GRB 050502B. A steep decay at the end of the plateau, as in example 4d (Figure [T4"]) . 
has been observed in the X-ray light curve of a few GRB afterglows such as GRB 060413 
and GRB 100508A. A fast decay at the end of the plateau, followed by a rebrightening, as 
in examples 6c and 6d (Figure [23]). has been observed in the X-ray light curve of a few GRB 
afterglows such as GRB 100814A, GRB 100418A, and GRB 060607A. 

This gives rise to interesting prospects of interpreting puzzling afterglow features within 
the RS scenario. In general, the observed afterglow emission should be a superposition of 
the FS and RS emission components. If the shock parameters of both shocks are the same, 
the RS component is usually over-shone by the FS component. However, since the FS is 
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ultra-relativistic while the RS is mildly relativistic, it is possible that the two shocks have 
different values of e e and/or e#. Moreover, as the ejecta may carry a primordial magnetic 
field from the central engine, it is reasonable to assume that the RS is more magnetized than 
the FS. If the magnetization of the ejecta is not too high (e.g. the Poynting-flux-to-matter- 
flux ratio o < 0.1), the magnetization of the ejecta would not affect the strength of the RS 
but would enhance synchrotron emission from the RS^|. As a result, the strengthened RS 
emission would become comparable or even brighter than the FS emission. So it is possible 
that the observed afterglow emission indeed includes contributions from both the FS and 
RS. 



6. Conclusion 

We have investigated in detail the dynamics and GRB afterglow light curves for a 
relativistic blast wave with a long-lived RS. This long-lived RS is formed in a stratified 
outflow ejected from a central source, which spreads out radially and forms various density 
structures in the flow due to its stratification. Due to spreading of Lorentz factor, this 
ejecta flow gradually catches up with the blast wave and adds its kinetic energy into the 
blast, naturally maintaining a long-lived RS wave. As a result, this blast wave is not in 
the self-similar stage, as is described in the BM76 solution. Instead, the blast wave is being 
continuously pushed by the RS, which displays various forms of energy injection scenarios. 

In order to find such dynamics of a blast wave with a long-lived RS, we make use of 
Ull with the mechanical model and perform detailed numerical calculations. Investigating 
a total of 20 different shapes of ejecta stratifications, we explain the effects and consequences 
of radial spreadings on the FS and RS dynamics. In particular, we show that there exists 
a whole new class of the RS dynamics with fast and strong evolutions. The FS dynamics 
is also shown to exhibit consistent behaviors for those diverse types of energy injections 
(Figure 125]) . A high accuracy shown in Figure |25] indicates that we are presenting here a 
"precision dynamics" for the blast waves with a long-lived RS. 

Employing a Lagrangian description of the blast wave, we perform a sophisticated calcu- 
lation of afterglows. In particular, our calculation has (1) a spatial resolution into the blast 
wave region and (2) a pressure profile that smoothly varies over the blast. For every shell 
on the blast, we keep track of an evolution of (1) the thermodynamic quantities of shocked 



8 In the regime of a higher a, the blast wave dynamics would be affected since the RS shock jump condition 
is significantly modified (Zhang & Kobayashi 2005). The calculations in this paper are relevant when such 
an effect is not important. 
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gas (pressure, energy density, adiabatic index, etc) and (2) the magnetic field and power-law 
spectrum of electrons. The FS and RS light curves are found by integrating over the entire 
FS and RS shocked regions, respectively, while making use of an analytic expression for 
observed spectral flux 5F° h * , which we derive here in terms of an observed frequency z/ t, s 
and observer time t D b s (Section l3.3p . 

The resulting afterglow light curves display interesting features. Since the FS strength 
mainly depends on the Lorentz factor of the blast, the FS light curves do not sensitively 
depend on the ejecta stratification. The strength of the RS, on the other hand, sensitively 
depends on the ejecta density flow p C j(RS) entering the RS, so that rich afterglow features 
show up in the RS light curves. As demonstrated with our 20 examples, the RS light curves 
naturally produce diverse and distinctive features (Figure I2"7j) as we move through the space of 
ejecta stratifications (Figure |2"41) . In particular, designing proper stratifications in the ejecta, 
the RS light curves reproduce many observed X-ray features, including various temporal 
breaks (Figures 151 and [IT]) and decay indices (Figures [T71 and 120]) . plateaus (Figures IB] [TT] 
and mj), steep declines (Figures [8] HU and [23]), bumps (Figure [TT]) . and re-brightenings 
(Figure [23]). Since the FS and RS could have different efficiency in particle acceleration and 
the GRB ejecta is likely more magnetized than the ambient medium, it is plausible that the 
RS emission would become as bright as or even brighter than the FS emission. Therefore, we 
believe that the RS could be a strong candidate to account for the observed GRB afterglows. 

ZLU is grateful to Massimiliano De Pasquale for helpful discussions. This work was 
supported by NASA NNX10AD48G, NSF AST-0908362, "Research in Paris" of the Paris 
City Hall, and CRI of MEST/NRF. 
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Fig. 1. — Illustrative diagram of a spherical blast wave. A forward shock (FS) wave sweeps 
up the surrounding ambient medium, and a reverse shock (RS) wave propagates through 
the ejecta. The shocked ambient medium is separated from the shocked ejecta by a contact 
discontinuity (CD). A Lagrangian description is employed to track an adiabatic evolution of 
all shells on the blast between the FS and RS. 
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Fig. 2. — Schematic figure of a spherical shell at radius r at time t expanding with a Lorentz 
factor r. A photon emitted in the direction of the observer at time t by an electron at a 
polar angle 9 is received by the observer at an observer time t Q ^ s as given in Equation ( 1201) . 
A thin ring between 9 and 9 + 59 is considered in order to derive an analytical expression 
for observed spectral flux 5F° h * ; see Equation fl38|) . 
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Fig. 3. — Ejecta stratifications for examples la, lb, and lc. The ejecta Lorentz factors 
r e j(r) are shown as a function of the ejection time r. 
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Fig. 4. — Blast wave dynamics for examples la, lb, and lc. In the panel (a), the curves 
denoted by T show the Lorentz factor of the blast wave as a function of the radius r r of the 
RS, and the curves by r e j(RS) show the Lorentz factor of the ejecta shell that gets shocked 
by the RS when the RS is located at radius ry. In the panel (b), the curves denoted by pf 
show the pressure at the FS, and the curves by p r show the pressure at the RS. The panel 
(c) shows the relative Lorentz factor 743 across the RS wave, as is given by Equation (III]) . 
The panel (d) shows the density n e j(RS) = p e j(RS)/m p of the ejecta shell, which enters the 
RS at radius r r . 
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Fig. 5. — Afterglow light curves for examples la, lb, and lc. The panel (a) shows the FS 
emissions in X-ray (1 keV) and R band as a function of the observer time £ b s . The panel 
(b) shows the RS emissions in X-ray (1 keV) and R band. 
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Fig. 13. — Same as in Figure HI but for examples 4a, 4b, 4c, and 4d. The F curves are 
multiplied by a factor of 0.7, in order to avoid an overlap with the r e j(RS) curves. 
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Fig. 16. — Same as in Figure HI but for examples la, 5a, and 5b. 
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Fig. 19. — Same as in Figure HI but for examples 6a, 6b, and 4d. 
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Fig. 21. — Same as in Figure |3l but for examples 6c, 6d, and 6e. 
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Fig. 25. — Lorentz factor F of the blast wave for all 20 examples. 
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Fig. 26. — Lorentz factor T of the blast wave for examples la and 5b. The T curves named 
by la and 5b are taken from Figure [161 The F curves denoted by la c and 5b c are found for 
examples la and 5b, respectively, by making use of a customary pressure balance p r = pf 
across the blast wave. 
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